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Abstract: We present a detailed study of the production of a high transverse- momentum 
lepton pair at hadron colliders, which includes the exact 0{q) electroweak corrections 
properly matched with leading logarithmic effects due to multiple photon emission, as 
required by the experiments at the Fermilab Tevatron and the CERN LHC. Numerical 
results for the relevant observables of single Z-boson production at hadron colliders are 
presented. The impact of the radiative corrections is discussed in detail. The presence in 
the proton of a photon density is considered and the effects of the photon-induced partonic 
subprocesses are analyzed. The calculation has been implemented in the new version of 
the event generator HORACE, which is available for precision simulations of the neutral and 
charged current Drell-Yan processes. 
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1. Introduction 

At hadron colliders, such as the Fermilab Tevatron and the CERN LHC, the production 
of a high transverse-momentum lepton pair, known as Drell-Yan process Q, plays an 
important role: it allows, in the neutral current channel, to study the physics of the Z boson 
and, in particular, to determine the effective weak mixing angle from the measurement of 
the forward-backward asymmetry Q; in the charged current channel, a high precision 
determination of two fundamental parameters of the Standard Model, namely the mass 
and the decay width of the W boson, can be obtained Q; it provides, both in the neutral 
and charged current channel, stringent constraints on the density functions which describe 
the partonic content of the proton it can be used as a standard reference process 
and, therefore, as a luminosity monitor of the collider [^, |6|. Furthermore, it represents a 
background to the search for new heavy gauge bosons ||7|. 

The accuracy in the determination of the theoretical cross section has greatly increased 
over the years. The calculation of next-to-leading order (NLO) QCD corrections Q has 
been one of the first test grounds of perturbative QCD. Next-to-next-to- leading order 
(NNLO) QCD corrections to the total cross section have been computed in ref. but dif- 



ferential distributions with the same accuracy have been obtained only recently in ref. |10| 



The size of the NNLO QCD corrections and the improved stability of the results against 
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variations of the renormalization/factorization scales raises the question of the relevance 
of the 0{a) electroweak (EW) radiative corrections, which were computed, in the neutral 
current channel, first considering the gauge-invariant subset of QED corrections |11] and 
then the complete set of 0{a) EW corrections |12, 13 1. 



A realistic phenomenological study and the data analysis require the inclusion of the 
relevant radiative corrections and their implementation into Monte Carlo event genera- 
tors, in order to simulate all the experimental cuts and to allow, for instance, an accurate 
determination of the detector acceptances. The Drell-Yan processes are included in the 
standard QCD Parton Shower generators HERWIG Q and PYTHIA Recently, there has 
been important progress to improve the QCD radiation description to NLO, which has been 
implemented in the code MC@NLO Another important issue is a reliable description of 
the intrinsic transverse momentum of the gauge bosons, which can be obtained by resum- 
ming up to all orders contributions due to multiple soft-gluon radiation. The generator 
RESBOS [17], used for data analysis at Tevatron, includes these effects. 

If the inclusion of QCD radiation is mandatory for the simulation of any process at 
hadron colliders, one should not neglect the impact of EW corrections on the precision 
measurement of some Standard Model (SM) observables, like the Z boson mass and de- 
cay width, important for detector calibration. For instance, the generator ZGRAD [11, 12 1 
includes the exact 0{a) EW corrections, which have been shown to induce a shift on the 
value of ruz extracted from the Tevatron data of few hundreds of MeV , depending on the 
experimental set up, mostly due to final-state QED radiation. Some event generators can 
also account for multiple-photon radiation: in the published version of HORACE |jl8|, 1£ , ^ 
QED radiation in leading-log approximation was simulated by means of a QED Parton 
Shower ||2^; the standard tool PHOTOS [^2| can be adopted to describe QED radiation in 
the Z decay. 

Since the Drell-Yan events can be used, in principle, to determine the collider lu- 
minosity at a few per cent level, the theoretical cross section must be known with the 
same accuracy, requiring also the inclusion of 0{a) EW corrections. Furthermore, the 
0{a) EW contributions give large corrections to the high tails of the invariant mass and 
lepton transverse momentum distributions, because of the presence of large EW Sudakov 
logarithms [^], thus affecting significantly the SM normalization in the search for new 
heavy gauge bosons. 

Neutral and charged current Drell-Yan processes can be combined to obtain an inde- 
pendent way of measuring the W boson mass: in fact, it has been shown [22] that the 
ratios of W and Z observables, properly defined, is less sensitive to missing higher-order 
QCD eff'ects and to the pdf uncertainties; on the other hand, these ratios are sensitive to 
the value of the W mass and can be exploited for its precise determination. The impact 
of the EW corrections in this framework is relevant. Preliminary studies ]24] have shown 
that these corrections do not cancel in the ratio and could play an important role in the W 
mass determination according to this method. It is therefore very important to have a tool 
to perform a realistic simulation of the Drell-Yan processes, including 0{a) and multiple 
photon corrections, not only in the charged current channel [^], but also in the neutral 
current case. 
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The aim of this paper is to present a precision EW calculation of the neutral current 
Drell-Yan process, which includes the exact 0{a) EW matrix elements properly matched 
with leading logarithmic higher-order QED corrections in the Parton Shower approach. The 
matching of perturbative corrections with Parton Shower, which is a topic of great interest 
in modern QCD simulations |25], has already been illustrated in ref. [20| and is realized 
along the lines already presented in ref. The use of the pdf set MRST2004QED pTj , 

which describes the partonic content of the proton also in terms of a photon density, implies 
the calculation of photon-induced partonic processes. The latter contribute to the inclusive 
cross section a (^pp — > l~^l~{n^) + and are of the same perturbative order as the Born 
Drell-Yan tree-level reaction (subprocess 77 l^l^) or of the same perturbative order of 
the 0{a) corrections (subprocess 75 — > l'^l~q). The latter have been recently addressed to 
in ref. @. 

Several distributions of physical interest are analyzed, disentangling the effect of differ- 
ent classes of radiative corrections and discussing various sources of theoretical uncertainty. 
The calculation is implemented in the new version of the Monte Carlo event generator 
HORACE *, which combines, in a unique tool, the good features of the QED Parton Shower 
approach with the additional effects present in the exact 0{a) EW calculation. This task 
is non trivial from several technical points of view and faces all the conceptual problems 
of developing a NLO event generator. 

The paper is organized as follows. In Section ^ we present the calculation of the 
partonic subprocesses and of 0{a) radiative corrections. In Section 3A we describe the 
matching of the fixed-order results with the QED Parton Shower. In Section |3^ we present 
the computation of the hadron-level cross section a (^pp l~^l~{n^) + X j and discuss 
the subtraction of the initial-state collinear singularities to all orders. In Section ^ we 
present phenomenological results for several physical distributions and discuss the impact 
of 0{a) EW corrections, photon-induced processes and higher-order QED contributions. 
Finally, in Section ^ we draw our conclusions and discuss possible developments of this 
work. 



2. Partonic process: matrix elements calculation 

We consider a description of the proton which includes also the presence of a photon as a 
parton, described by the corresponding density function, as done for instance in ref. p7| . 
The lowest-order cross section for the production of a high transverse-momentum lepton 
pair starts at 0{a'^) and receives contributions from two partonic subprocesses, namely 
q{pi) q{p2) ^"(Pa) ^^(Pi) and also 7(^1) 7(^2) — > ^"(Pa) ^^(^4)- When considering this fi- 
nal state at O(a^), we must include, in addition to the usual radiative process q{pi) q{p2) 
^ (Pa) ^^(^4)7(^5)1 also a new subprocess, i.e. q{pi) 7(^2) ^ (Pa) ^^(^4)9(^5)) which con- 
tributes to the inclusive signature. 

2.1 Born approximation 

We start considering the neutral current Drell-Yan partonic process q{pi) q{p2) ^ (Pa) ^^(P4); 
*The code can be downloaded from the url http://www.pv.infn.it/hepcomplex/horace.html 
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Figure 1: Born diagrams for the qq (a) and for the 77 (b,c) subprocesses. 



which is depicted in figure |l|(a). This process is a neutral current process and its ampUtude, 
neglecting the Higgs-boson contribution, is mediated by s-channel photon and Z-boson ex- 
change. In the unitary gauge, the tree-level amplitude reads as 



Mo = M^ + Mz (2.1) 



s 



2 ^ ^ g^^ k^k^/s 



S 

o2 



dfiu k^kp/ s 

z 



Mz = 2~2 1 I -n ["(^'2) {Vq f + 0,7^75) u{pi)\ [^(^3) {vi Y + 15) v{Pi)] 

sici s — mj, + iL zfUz 



° 5/^1^ kfikj^/ s 



sjcj s-ml + iTzTUz ^'^e 

where is the Z-boson mass and Tz is the Z decay width, necessary to describe the Z 
resonance region, s = {pi +^2)^ is the squared partonic center-of-mass (cm.) energy and 
k'^ = Pi +^2' ~ e^/(47r) is the fine structure constant, ce = mw/mz is the cosine of 
the weak mixing angle. The vector and axial-vector couplings of the Z-boson to fermions 
are Vf = Tj — 2QjSg and aj = —Tf where Tf = ±1/2 is the third component of the weak 
isospin and is the electric charge of the fermion /. 

The subprocess j{pi) 7(^2) ^"(Pa) ^^(^4)1 which is depicted in figure |l| (b,c), is, at 
lowest order, a pure QED reaction, whose differential cross section, in the partonic cm. 
frame and neglecting all fermion masses, reads as 

da.yy 2txo? f ^ + cos^ 6^^ 

— (2.2) 



d cos 9 s \ sin^ 
2.2 The 0{a) calculation 

The complete 0{a) EW corrections to the neutral current Drell-Yan process have already 
been computed in refs. |12, |^]. We have repeated independently the calculation and 



included in addition the photon-induced processes. We summarize here the main features 
of our approach. 

The C(a) corrections include the contribution of real and virtual corrections. The 
virtual corrections follow from the perturbative expansion of the 2^2 scattering amplitude 
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Figure 2: Some examples of one-loop virtual diagrams. 



M=Mo + A^r* + • • • and contribute, at 0{a), with 2Re{Ml'''^ M*^) . The 0{a) virtual 
amplitude includes two contributions, namely the one- loop renormalization of the tree- 
level amplitude and the virtual one-loop diagrams. The real corrections are due to the 
emission of one extra real photon and represent the lowest order of the radiative process 
q{pi)q{p2) — > {Pz)l^ {Pi)l{k) ■ They can be further divided into soft and hard corrections. 



Ml = M\°^* + Ml'''''. The for mer satisfies, by definition, the Born-like 2^2 kinematics 
and can be factorized as = ^sbIA^oPj where 5sB is a universal factor that depends 

only on the properties of the external particles. The total cross section includes soft and 
hard corrections and is independent of the cut-off used to define the two energy regions. 
Virtual and real soft corrections are separately divergent due to the emission of soft photons, 
but the divergence cancels in the sum of the two contributions. 

2.2.1 Virtual corrections 

The 0{a) virtual corrections to a 2 — > 2 reaction include the contribution of counterterms, 
self-energy, vertex and box corrections. Few diagrams representative of the different kinds 
of corrections are depicted in figure ^. The 0{a) virtual corrections have been calculated 
using the packages FeynArts and FormCalc p9|, The numerical evaluation of the 1-loop 



integrals has been done using the package LoopTools2 [30|, based on the library ff pl| . 
We will write the 1-loop virtual amplitude as M^""* = M^J' + Aif^ + TW^*^*^*^^ + 
where M^^ includes all the counterterms and the wave function corrections on the external 
legs, M'S^^ describes the self-energy corrections to the photon and to the Z propagator and 
the contribution due to the ^ — Z mixing and M^'^^^^' °^ describe respectively the vertex 
and the box corrections. The mass of the fermions in the scalar 1-loop integrals regularizes 
in a natural way the mass singularities due to the emission of a (virtual) collinear photon. 
The infrared divergence of the integrals has been regularized by means of a small photon 
mass A. 
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The introduction of the Z decay width in the propagator of the Z boson is mandatory 
to describe the resonance region and to regularize the divergence due to the pole of the 
propagator. In the 't Hooft-Feynman gauge, the 1-loop corrections to the Z propagator 
read 

{-ign r^^T^ (Re(nzz(s)) + bm\ + (s - -^-^ (2.3) 

s — mi, + il zUiz s — rfiy 



.2 



where Ilzz{s) is the transverse part of the Z self-energy corrections, drriz and 6Zz are the 
Z mass and wave function renormalization constants, respectively. The two counterterms 
cancel the divergences present in the self-energy corrections. We remark that we have to 
take the real part of Iizz{s) in order to avoid a double counting with the decay width in the 
denominator. In fact the 0{a) contribution is 2Ke{Ai'^^^J Mq): the real part prescription 
in the latter expression avoids the double counting in the interference with the Z tree- 
level diagram, but not the one in the interference with the photon tree-level diagram. 
Furthermore, the second factor — m|) in eq. (|2.3|) is not corrected by the decay width, 
again to avoid a double counting; we can check, by expanding the self-energy corrections 
around s = m^, that the 0{a) expression is regular for s — > m|. 

The 'J — Z mixing diagrams require the introduction in the Z propagator of a decay 
width to describe the Z-resonance region, i.e. 

Hg^n fRe(n,^(.)) + i^^^6Z,z + bzz,) - (2.4) 

s — niz + il z^z \ ^ 1 J s 

Also in this case we need to take the real part of 11^^ to avoid the double counting between 
the imaginary part of the self-energy and the imaginary term proportional to the decay 
width. 

Examples of vertex and box corrections are depicted in figure |2[ The abelian vertex 
diagrams and the box diagrams with a photonic correction are infrared divergent. All 
the vertex and box diagrams with a Z-boson exchange yield the so-called EW Sudakov 
logarithms, namely terms like alog^ (s/m|), whose importance grows for large invariant 
mass of the final-state lepton pair, while they are almost negligible at the Z resonance. The 
box diagrams with a photon and a Z boson in the loop yield a logarithmic divergence at 
s = m?z- Since all the terms with log(s — m|) form a gauge invariant subset of corrections, 
we regularized the divergence by replacing everywhere log(s — m|) — > log(s — ml + ir^m^). 
The 7Z boxes are resonant at s ~ m|. The divergence has been regularized by modifying 
the box amplitude in the following way 

2 

^ Pt"^ + ^-^^ A^'f + A^ll (2.5) 

' ' s — rriz + li- zTTLz 

The calculation has been repeated with two different gauge choices, namely the 
with ^ = 1 gauge and the background field gauge, with parameter g = 1 |3^ ]. The two 
results perfectly agree, and this is an important check of the calculation of the bosonic 
self-energy and of the non-abelian vertex corrections. 
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Concerning the renormalization of the 1-loop amphtude, the ultraviolet divergences 
which appear from the virtual diagrams can be cancelled with the mass, Sm'^, and wave 
function, 5Zz, renormalization constants of the Z boson, with the wave function renormal- 
ization of the photon, 6Z^^, and of the jZ mixing, 5Z^z and 5Zz'y and by the renormal- 
ization of the two vertices Zqq and Zl^l^ . The latter include the charge renormalization 
and the wave function renormalization of the external fermions, of the photon and of the 
Z boson. We use the electric charge and the gauge boson masses e,mw,'n^z as input pa- 
rameters, we replace the bare coupling cq in terms of the renormalized quantity and of the 
corresponding counterterm cq = e{l — 5e/e). The electric charge counterterm is fixed by the 
request that in Thomson scattering the renormalized charge is given by the fine structure 
constant; its expression depends on the quark masses running in the photon vacuum po- 
larization, the value of which can be adjusted in order to make the running electric charge 
reproduces the value a(m|) ||3^. For the sake of brevity, we can introduce a counterterm 
for the weak mixing angle, which is given by a combination of the mass counterterms of the 
W and Z bosons, following from its definition. The W- and Z-boson mass and wave func- 
tion renormalization constants, the 77 and 7Z wave function renormalization constants 
are defined in the on-shell scheme. 

The choice of the input parameters of the SM lagrangian has a numerical impact on 
the calculation of the physical observables. In order to obtain predictions with the smallest 
possible parametric uncertainty, it would be convenient to choose as input parameters 
a, G^,mz, which was the common choice at LEPl. The drawback is that in this approach 
m^, is a prediction of the SM, whereas, at hadron colliders, it would be interesting to keep 
it as an input, which can be measured fitting the data. It is therefore more useful to choose 



a,mw,'mz as inputs. For instance, the method proposed in ref. [23| to study the ratio of 
W and Z observables and to use it to determine the W mass requires m^i, as input of both 
charged and neutral current simulations, which will then be compared with the data. 

Unfortunately, the choice a, , suffers, if strictly applied, of the ambiguity due to 
the choice of the quark mass values in the photon vacuum polarization. The ambiguity can 
be removed by re-expressing, in the Born amplitude, the fine structure constant in terms 
of a different quantity, like the Fermi constant or the effective electromagnetic coupling 
at the scale g^, a(g^), whose 0{a) expressions exactly cancel the quark mass dependence 
due to the vacuum polarization diagrams. In order to rewrite the tree- level amplitude in 
terms of these new couplings, we follow the prescriptions developed at LEPl to study the 
Z-resonance and treat separately the photon- and the Z-diagram contributions. In we 
replace 

- eHq') = / (l - Aa{q')) (2.6) 
Aa(g2) = n(/)(g2)-n(/^)(0) 

where Aa{q'^) expresses the running of the electric charge due to the fermionic photon 
vacuum polarization (nl*(?(g^) is related to the transverse part of the photon self-energy 
by Ayiy {q'^) = q'^IVyil {q'^) ). The quantity Aa includes also the non-perturbative hadronic 
contribution due to a loop of quarks at low virtualities. For the evaluation of Aa^^ we 
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(a) (b) (c) (d) 

Figure 3: 0{a) bremsstrahlung diagrams. 



use the routine hadrBn |^^. Because we include the photon vacuum polarization effects 
in the lowest-order coupling, we have to subtract the 0{a) expansion of e^(g^), to avoid a 
double counting when we include the full set of 0{a) corrections. 

In the case of M.z we can rewrite e^/ (sgC^) as ff^/cg and then use the relation, computed 
up to 0(0), of the weak coupling g with the Fermi constant and the VF-boson mass 

The quantity Ar represents all the radiative corrections to the muon-decay amplitude 
2.2.2 Bremsstrahlung corrections 

The real radiative corrections to the neutral current Drell-Yan process, described by the 
amplitude A4 1 , are given by all the Feynman diagrams (figure ^) with the emission of one 
extra photon from all the electrically charged legs of the Born diagrams. 

The probability amplitude has been calculated in the unitary gauge with massive 
fermions. We integrate the squared matrix element over the whole photon phase space 
and split the allowed photon energy range in two intervals, [A, AS] and [AE, E^ax]- The 
cut-off AE <C -y/s is chosen in such a way that the photon with smaller energy is considered 
soft and does not modify the 2 — > 2 kinematics of the Born amplitude. The small photon 
mass A has been introduced to regularize the infrared divergence. In this energy region the 
phase space integral, including the full angular integration, can be solved analytically. The 
result can be expressed in a factorized form, as 



{2-k)^2E. 



^ \Mi\^ = \Mo\' V 5sB{f,X) (2. 



f=q,q,l+l- 



where the soft Bremsstrahlung factor, see e.g. ref. [^], depends on the mass and electric 
charge of the external radiating particles and the phase-space region is defined by the 
request that the photon energy E^ satisfies \ < E^ < AE. We have explicitly checked 
that the sum of the virtual and soft-real contributions is independent of the choice of the 
photon mass A, in the limit of small A values. 

In the hard energy region the phase-space integration has been performed numerically, 
with Monte Carlo techniques improved by importance sampling to take care of collinear 
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Figure 4: Photon-induced process diagrams. 



and infrared singularities, as well as the peaking behaviour around the Z resonance. The 
sum of the soft and of the hard photon cross sections is independent of the cut-off Ai?. 
We have checked the independence of our numerical results from the choice of the infrared 
separator e = /\E/E for 10"^ < e < 10"^. 



2.2.3 Photon-induced processes 

In ref. |27| it has been proposed a new parametrization of the partonic content of the 
proton, which also includes a photon probability density. When using this set of pdf, the 
inclusive cross section a (^pp l^l^ + receives contributions also from the partonic 

subprocesses ^(^1)7(^2) {Pz)l~ {P'i)<l{k) (photon-induced), depicted in figure The 

latter are of the same perturbative order as the real bremsstrahlung corrections described 
in the previous subsection, i.e. they are an 0{a) correction to the Born process of eq. (p.l|). 
The squared amplitude of the photon-induced processes can be obtained by crossing sym- 
metry from the real bremsstrahlung one, evaluating the latter with the exchange {p2 <^ —k) 
and multiplying the result by a (—1) factor to account for the exchange of a fermionic line. 



2.3 Higher-order electroweak effects 

To incorporate higher-order EW corrections in a Born-like expression written with effective 
couplings, we followed the approach of ref. |3^, where the tree-level amplitude has been 
improved and takes into account all the self-energy and vertex corrections. The latter have 
been included by defining an effective overall coupling and an effective weak mixing angle. 
The amplitude Aiz becomes 

^ i8 G^ml pfi[q^) Jz,qq ■ Jz,i+i- ^2 9) 

^/2 1 - 5pirr - m\ + iTzUlz 

where the coupling vj of eq. ( |2.1D is replaced hy Vf = Tj — 2Q fKf{q^)s1. The definition 
of the quantities pfi, 6pirr, Kf{q'^) can be found in ref. ||3^. Eq. ( |2.9| ) incorporates also 
higher-order effects beyond 0{a), because of the resummation of 5pirr and of the fermionic 
part of the Z self-energy contained in pfi. Furthermore, 6pirr = SpUl + ^Pirl contains also 
leading two-loop corrections. In the amplitude A^-y we replace the fine structure constant 
with the running electromagnetic coupling according to eq. ( |2.6| ). 
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For the numerical analysis we combine the exact 0{a) corrections described in the 
previous subsections and evaluated with a{0),mw,'mz as input parameters, with the Born- 
like expressions of this section. In order to avoid a double counting in the combination, we 
subtract the 0{a) content of eq. ( |2.9| ). All the numerical results shown in the following are 
obtained according to the EW scheme here described. 



3. Matching QED higher orders and hadron-level cross section 
3.1 Matching 

In this section we describe the matching of the 0{a) EW calculation with higher-order QED 
corrections (cfr. ref. (2^). The latter can be included in a generic scattering cross section 
in the QED Parton Shower approach, which resums to all orders the leading logarithmic 
effects. At 0{a) the Parton Shower reproduces only the QED leading- log approximation 
of the exact 0{a) EW calculation presented in section |2[ We would like to combine the 
exact 0{a) results and QED higher orders, to improve the approximation intrinsic to the 
Parton Shower, avoiding at the same time double counting at 0{a). A detailed description 



of the matching procedure can be found in refs. 26|. Here we recall the basic ideas and 
results. 

Our master formula is 



= F5v n(g2,e) ^- n^^^'O \^n,LL\'' d<!>n (3.1) 

n=0 ' \i=0 / 

and can be understood in the following way: 1) the tree-level cross-section (n = in 
the sum) is corrected, in leading-log approximation, to all orders by the Sudakov form 
factor n(Q^,e), which accounts for virtual and soft-photon emission up to a scale e in a 
hard-scattering process characterized by a virtuality scale Q; 2) the resulting expression is 
dressed by the QED Parton Shower, with the real bremsstrahlung n-photons squared matrix 
elements in leading-log approximation; 3) the correction factors Fsv and Fn^i provide the 
remaining 0{a) corrections missing in the leading-log approximation: in particular, Fsv 
contains the remainder of the soft plus virtual corrections, whereas Ffj^i gives, for the 
real photon emission, the correction due to the exact bremsstrahlung matrix element with 
respect to the Parton Shower approximation. 



The expansion at 0{a) of eq. (3T) coincides with the exact NLO cross section. Fur- 
thermore, all higher-order leading-log contributions are the same as in a pure QED Parton 
Shower. It is worth noticing that Fsv^Fn^i are, by construction, infrared safe and free of 
collinear logarithms. 

3.2 Hadron-level cross section 

In this section we discuss the calculation of the hadron-level cross section a{pp — > + 
nj + X) and the procedure to subtract the initial-state mass singularities. This procedure 
makes the resulting hadron-level cross section independent of the (fictitious) value of the 
initial-state quark masses. 
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The subtraction at 0{a) has been discussed, for example, in ref. |37] and is obtained 
by a redefinition of the parton densities. The hadron-level cross section at 0{a) can be 
written as 



da{pp 



i+r + x) 



E 

a.,b=q,q,'r 



dxidx2 qaixi,M^)qh{x2,M'^ 



daf 



+ da; 



ab 



(3.2) 



- {/^qa{xi,M^)qdx2,M^)+qa{xi,M'^)Aq^{x2,M^)) da^ 



where a,b run over all parton species described by the densities qi{x,M'^), M is the fac- 
torization scale, diTQ* and da^ are the Born and 0{a) partonic cross sections, initiated by 
partons a and b, respectively. In the case of photon-induced processes, da^'' = da^"^ = ^. 
The relevant 0{a) subtraction terms are 



Aq,{x,M^ 



dz qil-,M']—Q. 



a 



2tt 



Pr 



+ fq{z) + q-y 



.f2 



a 



(-) ( log ( ^ 



27r 



Qi 



p. 



(^) lo: 



21og(l - z) - 1 



M2 
2 



m 



(3.3) 



Aq^{x,M'^ 



E 

i=q,q 



dzq,(-,M']—Qi 



a 



Pr 



q-*iq\ 



(z) log ^ -21og(l-z)-l 



where Qi and rrij are the electric charge fraction and the mass of the quark i\ the functions 
fi{z) {i = q,'y) allow to change the subtraction scheme (e.g. DIS or AIS) and are defined 
as [H 



fqiz) 



log(l - z) 
1 - z 

+ 3 + 2z- 



1 



9 vr 
2 ^ "3 



2(l-z)_ 

2^ 



:i + 2)iog(i 



l + z^ 



log(z) 



f^{z) = {{I -zf + z^) log 



1 



iz^ + 8z-l 



(3.4) 



Since also the photon-induced processes contribute to the hadron-level cross section and 
develope a mass singularity when the outgoing quark is collinear to the incoming photon, 
also the Altarelli-Parisi splitting function P-y-^qq contributes to the subtraction term for the 
quark densities. The processes jq l^l^q, because of a photon exchange in the t-channel 
in the peripheral diagram, develop a collinear singularity, proportional to the splitting 
function Pq^^q(z), which can be reabsorbed in the photon density. In the DIS scheme, the 
function f{z) can be fixed by any combination of observables and will be therefore set to 
zero. 

Given the presence in the hadron-level cross section of eq. (|3.2| ) of the product of two 
parton densities, the subtraction procedure in a factorized form could yield terms of 0{a^), 
which have been discarded at 0{a) for consistency. 



^In the numerical evaluation, we omit 0{a) corrections to the subprocess da'^'^ , because numerically 
negligible. 
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The generalization of the independence from the value of the quark masses of the 
0{a) cross section of eq. ( p.2[ ) to the cross section including also QED higher-order correc- 



tions has been discussed in detail in ref. |2C] and the main result is our master formula for 
the computation of the hadron-level cross section and event simulation 



dcThad = / dxidx2 qa{xi,M'^)qb{x2,M'^) X 



a,b=q,q 



I FsvUiQ^e) (IIFhA \Mn,LL\ 

[ n=0 \i=0 / 



d<S>r 



+ 



d&. 



ab,: 



Aqa{xi,M'^) , Aqb{x2,M 



+ 



qaixuM^) (?b(x2,M2) 



da. 



ab 



where daU^j^jdaJl^ can be derived straightforwardly from eq. (| 

ab,sub 



(3.5) 



(3.6) 



The factors with a tilde and also daa ' represent quantities subtracted of the initial- 



state singularities (cfr. Section 4 of ref. |2C] for the definitions and more details) 



4. Numerical results 



All the numerical results have been obtained using the following values for the input pa- 
rameters 



a = 


1/137.03599911 




= 1.16637 10-5 GeV-2 


rriz = 


91.1876 GeV 


rriw 


= 80.398 GeV 




= 2.4952 GeV 


niH = 


115 GeV 


nie 


= 510.99892 KeV 




= 105.658369 MeV 


nir = 


1.77699 GeV 


rriu 


= 66 MeV 


rUc 


= 1.2 GeV 


nit = 


170.9 GeV 


rrid 


= 66 MeV 


rus 


= 150 MeV 


mb = 


4.3 MeV 


Vud 


= 0.975 


Vus 


= 0.222 


Vub = 





Vcd 


= 0.222 


Vcs 


= 0.975 


Vcb = 





Vtd 


= 


Vts 


= 


ytb = 


1 



and have been computed in the EW scheme defined in Section The set of parton density 
functions used to compute all the hadron-level cross sections is MRST2004QED In this 

set of pd/the QCD and the QED factorization scales are set to be equal and, as usually done 
in the literature [ll|, 12], we choose M = niz- The use of the pdf set MRST2004QED implies 
that our numerical results correspond to the DIS factorization scheme. The computation of 
the hadron-level results requires the numerical evaluation of the subtraction term defined 
in eq. (|3.3|); a grid of values in the variable x, which is then interpolated, is obtained by 
means of the numerical library CUBA [39|. All the hadron-level results refer to the LHC, at 
a nominal cm. energy ^/s = 14 TeV. 

The following cuts have been imposed to select the events 



P±,e > 25 GeV, |%| < 2.5 



(4.1) 
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Figure 5: Invariant mass distribution around the Z peak (left) and relative effect of different 
contributions (right), for bare muons and reconibined electrons. 




le-l)9 \ , , , , , \ I 1 , , , , 
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Jl/,+,- (GeV) J/,*,- (GeV) 



Figure 6: High tail of the invariant mass distribution (left) and relative effect of different contri- 
butions (right), for bare muons and recombined electrons. 

where p^^i and r]i are the transverse momentum and the pseudo-rapidity of the charged 
leptons, respectively. 

Our results are obtained for muon pair final states. However, we also show results for 
recombined electrons in the case of the invariant mass distribution. In fact, we assume 
perfect isolation of photons from the muon, which is experimentally achievable with good 
accuracy: the resulting correction is therefore amplified by large muon mass collinear 
logarithms, because the photon emission is not treated inclusively in the region around the 
muon. In the case of electrons, it is not possible experimentally to separate them from 
the photon track, when the latter lies within a cone around the lepton smaller than the 
detector angular resolution. We adopt the following recombination procedure 

• photons with a rapidity jry^j > 2.5 are never recombined to the electron; 

• if the photon rapidity is \r]^\ < 2.5 and Re-y = ^ {rje — rj^Y + (p^-y < 0-1 i4'e-y is the 
angle between the photon and the electron in the transverse plane), then the photon 
is recombined with the electron, i.e. the momenta of the two particles are added and 
associated with the momentum of the electron; 
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1. lowest order (Born), as in Section |2.3| , only qq subprocess 

2. exact 0{a) EW corrections to the qq subprocess 

3. exact 0{a) EW corrections to the qq subprocess plus 77 and ^7 contributions 

4. exact 0{a) EW matched with higher-order QED corrections to the qq subprocess 

Table 1: Different approximations for the calculation of the neutral current Drell-Yan cross section. 





1. 


2. 


3. 


4. 


a(pb) 


739.1(2) 


710.7(1) 


715.8(1) 


712.8(2) 



Table 2: Cross sections, in pb, using approximations 1., 2., 3. and 4. given in table 1 and according 
to the EW scheme defined in section 2.3. 




Jl/x (GeV) A/x (GeV) 



Figure 7: Transverse mass distribution around the Z peak (left) and relative effect of different 
contributions (right). 

• the resulting momenta should satisfy the cuts of eq. ( |4.1| ). 

In order to study the different effects of the radiative corrections and partonic sub- 
processes on the relevant observables, we will distinguish the approximations described in 
table |. 

In table § we compare the values of the total cross section, within the cuts specified 
above, calculated in the approximations 1., 2., 3. and 4. of table |l[ The effect of the 
0{a) corrections is negative, of about 4% of the lowest-order result, while photonic sub- 
processes and QED higher orders enhance the 0{a) cross section of about 1% and a few 
per mille, respectively. 

In figure ^ we show the invariant mass distribution of the final-state lepton pair, which 
is peaked at the Z-boson mass value and which, taking as a reference the high-precision 
LEP measurement, can be used to calibrate the LHC detectors. The radiative corrections 
significantly modify the shape of this distribution, as can be seen from the right panel of 
figure showing the relative effect, expressed in units of the Born cross section (approx- 
imation 1.), of the different contributions corresponding to approximations 2., 3. and 4. 
of table 1. There is a well-known effect p| , |T9| due to final-state photon radiation which 
increases, in the muon pair final state, by almost 80% the lowest-order result in the region 
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Figure 10: High tail of the lepton transverse momentum distribution (left) and relative effect of 
different contributions (right). 



below the Z peak. Around the peak, instead, the correction is negative of about —15%. 
The size of these corrections is significantly reduced when considering electron final states, 
in agreement with the results known in the literature |11] and due to the photon recom- 
bination procedure, which implies a partial cancelation of the mass logarithms, according 
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Figure 12: Pseudo-rapidity distribution of the final state lepton (left) and relative effect of different 
contributions (right). 



(1.1 




Figure 13: Forward-Backward asymmetry as a function of the invariant mass of the lepton pair 
according to the approximations 1., 2. and 3. of table 1. 

to the KLN theorem. The impact of the photon-induced processes is of the order of a few 
per cent and is particularly evident away from the resonance region. The effect of multiple 
photon emission amounts to a few per cent (see the inset of the figure) and can not be 
neglected for an accurate detector calibration and precision Z-physics studies at hadron 
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colliders, as already remarked in ref. [19|. The Z-boson mass shift induced by QED higher 
order corrections in fits to the invariant mass distribution has been already quantified in 
ref. and found to be about 10% and of opposite sign of that caused by one-photon 
emission. 

In figure ^ we show the invariant mass distribution in the high tail region, where the 
Drell-Yan processes represent an important background to the searches for new heavy gauge 
bosons. The distribution receives large negative radiative corrections, of the order of a few 
ten per cent, due to EW Sudakov logarithm, as previously emphasized in refs. [12, 13| ^. 
With electron final states, the photon recombination procedure reduces the negative effect 
of the 0{a) corrections. The multiple photon emission yields a few per cent effect in 
this invariant mass region, while the photon-induced processes raise the cross section by 
approximately +12% of the Born approximation. The 77 subprocess contributes 

with about -|-3% of the qq Born cross section. Actually, its contribution is suppressed by 
the smallness of the photon density and can be mainly observed at large invariant masses 
of the lepton pair. 

In figure we consider the distribution of the Z transverse mass, defined as 



M± = ^2p^^i+ p^^i- (1 - cos 4>i+i-) (4.2) 

where (pi+i- is the angle between the leptons in the transverse plane. We remark the 
different size of the effect of the 0{a) corrections, with respect to the invariant mass 
distribution, to this quantity, that can be useful, in association with the W transverse 
mass distribution, to measure the M^-boson mass. In the vicinity of the Z resonance, the 
relative contribution of 0{a) EW corrections amounts to ~ —17% and is about a factor 
of two of the corresponding effect on the W transverse mass distribution in the charged 
current channel [pO| , ^ . This can be simply understood in terms of the dominance of final- 
state QED radiation within the full set of 0{a) EW corrections around the Z resonance 
and, consequently, because of the presence of two radiating leptons in the neutral current 
channel. Photon- induced processes reduce the relative size of 0{a) EW corrections and 
give a contribution of some per cent, especially above the Z peak, while multiple photon 
corrections contribute at some per mille level. 

In figure |8| we show the high tail of the transverse mass distribution, which also rep- 
resents an important background observable to the searches for new neutral gauge bosons. 
The presence of the EW Sudakov logarithms reduces the cross section by 18 to 35% for 
1 TeV < M± < 4 TeV. In contrast to the invariant mass case, the contribution of the 
photon induced processes is smaller in this case, of the order of +7%, while the effect of 
higher-order QED corrections is of the same size as that observed in the high tail of the 
invariant mass distribution. 

In figure ^ the lepton transverse momentum distribution is shown around the Z peak 
and, in figure in the high momentum tail. The effect of the radiative corrections, as 



■'■It has been recently noticed that -weak boson emission diagrams, that contribute at the same order 
of 0{a) EW corrections, lead to a partial cancellation of EW Sudakov corrections, when the weak boson 
decays into unobserved vu or jet pairs ||40| . 
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well as of photon-induced processes, shows a pattern similar to the one observed for the 
transverse mass distribution and discussed above. 



The rapidity distribution of the Z boson is presented in figure 11, The 0{ol) EW 
corrections are negative, of the order of -4%, and almost constant in the interval \yz\ ^ 1, 
and are smaller in size, at the 2% level, for larger values of the rapidity. It is worth noting 
that 0{oi) EW contributions to such an observables are of the same order of magnitude 
as NNLO QCD corrections |10| and, therefore, both the effects need to be taken into 



account in precision measurements of the Z rapidity distribution. We notice the overall 
positive correction due the photon-induced processes, at the 1% level. The multiple photon 
emission has an almost negligible impact on this observable, at the per mille level, as already 



discussed in ref. |19 |. 



The numerical results for the pseudo-rapidity of the final state lepton and the relative 
effect of the radiative corrections and photon-induced processes on this distribution are 



illustrated in figure 12, showing a pattern quite similar to that observed for the Z rapidity. 

At hadron colliders it is possible to define a forward-backward asymmetry ApB and 
to derive from it a measurement of the leptonic effective weak mixing angle sin^ Q^^^ . The 
forward-backward asymmetry can be written as 



F{Mi+i-)= [ dcose* B{Mi+i-)=[ dcosr 

Jo dcosO* y_i 



dcosO* 



where 



cos ( 



= / , = [p+(r)p-(/+)-p+(/+>-(r)] (4.4) 

M{i+i^)^A'p{i+i-) + viiin-) 

= -^iE±p,), f = 1 (Tevatron), / = (LHC) (4.5) 

V2 Pzil^l ) 

and M(l^l^) is the invariant mass of the final-state lepton pair, p±{l^l^) and pz{l^l^) are 
the total transverse momentum and total longitudinal momentum of the pair, respec- 
tively. The asymmetry can be expressed in terms of sin^^^jj, with good approximation, 
as ApB = b(a — ^^^^^eff)- ^ detailed description of the effect of the 0{a) corrections 
on the coefficients a,b can be found in ref. |l^, together with an analysis of the relevant 
backgrounds. In figure ^ we present the asymmetry distribution, evaluated with the cuts 
of eq. ([4.1|), according to the approximations 1., 2. and 3. of table 1. It can be seen 
that 0{a) EW corrections are relevant and modify the shape ApB below the Z peak, 
whereas the photon-induced processes, as for the multiple photon corrections [|l^, do not 
contribute significantly to this observable. 

5. Conclusions 

In this paper we have presented a precision electroweak calculation of the neutral current 
Drell-Yan process. The theoretical approach is based on the matching of exact 0{a) EW 
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corrections with QED Parton Shower, to account for the effect of multiple photon emission. 
The use of the proton pd/ parametrization MRST2004QED allows for the presence of a photon 
density in the proton and gives rise to additional partonic subprocesses which contribute 
to the inclusive Drell-Yan signature. The impact of the different contributions studied 
in the paper turns out to be important, in association with higher-order QCD effects, in 
view of the precision measurements of EW parameters at the Fermilab Tevatron and the 
CERN LHC, as well as to validate the existing pd/ parametrizations and to assess the SM 
normalization in the search for new neutral heavy gauge bosons. 

The calculation has been carried out using as inputs a(0), m^^, rather than a LEP- 
like choice a(0), G^,mz, in order to keep both gauge boson masses as free parameters which 
can be fitted from the data, since W- and Z-boson physics are intimately related at hadron 
colliders. 

Possible perspectives of the present paper include tuned comparisons with independent 
calculations, in the spirit of the work done for single W production during the Les Houches 
2005 |41] and the TeV4LHC |42| workshops, a combination of EW and QCD corrections 
at the event generator level and the addition of higher-order contributions such as 2-loop 
EW Sudakov logarithms |43]. 
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